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Number-resolved imaging of excited-state atoms using a scanning autoionization microscope 
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We report on a scanning microscopy technique for atom-number-resolved imaging of excited-state 
atoms. A tightly focussed laser beam leads to local autoionization, and the resulting ions are counted elec- 
tronically. Scanning the beam across the cloud builds up an image of the density distribution of excited 
atoms, with access to the full counting statistics at each spatial sampling point and an overall detection 
efficiency of 21 %. We apply this technique to the measurement of a spatially inhomogeneous electric 
field with a spatial resolution of 50 /xm and a sensitivity to electric field gradients of 0.04 V cm^^. 

PACS numbers: 32.80.Ee, 32.80.Rm, 032.80.Zb 



Ultracold atomic gases constitute an almost ideal testbed 
for studying complex quantum phases of matter [1]. State- 
of-the-art experiments are now able to image the transition 
from a superfluid to an insulating state with single-atom 
resolution flUI]- Recently it has become clear that Rydberg 
states can be exploited to introduce tunable, long-range in- 
teractions, leading directly to the creation of many-body 
quantum states [4] that are associated with strong, long- 
range spatial correlations With control over the ge- 
ometry and excitation parameters, more exotic states are 
possible, where the Rydberg excitations form a crystal with 
true long-range order 

To probe these strongly correlated systems effectively, 
Rydberg atoms must be detected with single- atom sensi- 
tivity and micron-scale spatial resolution. Standard tech- 
niques, such as absorption and fluorescence imaging, have 
been a ppli ed to the detection of ions ifioll and Rydberg 
atoms 111 in in an ultracold plasma, but single-atom sen- 
sitivity is challenging as only a few photons are scat- 
tered by each atom. If the atoms are confined to individ- 
ual trapping sites, single Rydberg excitations can be de- 
tected via trap loss 11121 [l3l|. Such an approach was re- 
cently employed to measure the spatial correlation func- 
tion in a strongly interacting Rydberg gas [6]. Spatial 
correlations have also been observed using ion imaging 
techniques ISO, but stray electric fields from the ion op- 
tics can lead to unwanted Stark shifts that preclude the 
observation of crystalline states. Other proposals include 
the extension of electromagnetically-induced transparency 
techniques lfl4l [isil to detect single Rydberg excitations 

Mm. 

Here we report on a novel technique for the measurement 
of the spatial distribution of excited-state atoms, based 
upon core excitation ifisll of a two-electron atom. If one 
valence electron is in an excited state (typically principal 
quantum number n > 10), excitation of the second valence 
electron leads to rapid and highly efficient autoionization. 
Using a tightly focused autoionization laser, we selectively 
ionize the Rydberg atoms in a small spatial region. The 
resulting ions are counted using an electronic detector (see 
Fig. [U, with an overall efficiency of up to 20 %. By scan- 
ning the laser across the cloud an image of the spatial dis- 
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1 . (Color online) (a) Relevant energy levels and transitions 
(b) Overview of the experiment. The distribution of 
Rydberg atoms is defined by the tightly focused A2 beam, (c) 2D 
image of the Rydberg density distribution at n = 56. The y-axis 
projection shows a single slice (black points). The solid red line 
is the corresponding slice through the ground-state fluorescence 
image, with the amplitude scaled to give the best overlap. The x- 
axis projection shows a slice through the centre of the distribution 
overlaid with a Gaussian fit (red solid line). 



tribution of the Rydberg atoms is built up, with single-atom 
sensitivity and full counting statistics at each position. Cru- 
cially, this technique operates independently of any exter- 
nal fields, such as trap potentials or electric fields. We ap- 
ply this technique to the spatially resolved measurement of 
an inhomogeneous electric field, achieving a sensitivity to 
electric field gradients of 0.04 V cm~^ with a spatial reso- 
lution of 50 /xm. 

The extracted electric field compares well to other Ry- 
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dberg static electrometry techniques fi3\ which quote a 
sensitivity of ~ 100 mV cm~^ with a spatial resolution 
of 7 /xm. Recent experiments in strontium optical lattice 
clocks indicate that DC Stark shifts can lead to significant 
systematic errors |ll9|]. Mapping the electric field at the 
10 mV cm~^ level using Ry dberg states could reduce this 
uncertainty by a factor of ~ 100. We note that Rydberg 
excitation in the lattice clock has already been proposed as 
a sensitive black-body thermometer i2(M . 

To create the initial cold atom ensemble, ^^Sr atoms 
are cooled and trapped in a standard magneto-optical trap 
(MOT) operating on the 5s^ ^Sq — > 5s5p ^Pi transition at 
461 nm (Ai in Fig. [T](a)), loaded from a Zeeman-slowed 
atomic beam. A fluorescence image of the cloud is used 
to obtain the ground- state atom number and density distri- 
bution. Atom numbers of (5 ± 0.5) x 10^ are achieved 
at densities of (2 ± 0.2) x 10^ cm~'^ and temperatures of 
5 ± 0.5 mK. 

The cold atoms are excited to a 5snd Rydberg state 
usinga two-photon coherent population trapping (CPT, 
l 2lL I22II ) scheme (Ai and A2 in Fig. [T] (a)). The exci- 
tation lasers are counter-propagating to reduce the resid- 
ual Doppler broadening and are pulsed on simultaneously 
for 1 /xs. A 3 ± 0.2 G quantization field is applied, and 
both excitation beams are circularly polarized such that the 
mj = +2 Rydberg state is excited. The lasers are sta- 
bilized on resonance with their respective transitions us- 
ing sub-Doppler frequency modulation spectroscopy (Ai 
1 23 1), and electromagnetically-induced transparency (A2 
|2^|25|]) in a dispenser-based vapour cell [26]. 

The CPT spectrum is measured in the cold atoms by 
stepping the frequency of Ai, with A2 held on resonance. 
For the 5s56d ^D2 Rydberg state, we obtain a full-width 
at half maximum (FWHM) of 3.7 ± 0.2 MHz in the low- 
intensity limit, compared to the 5s^ ^Sq — > 5s5p ^Pi transi- 
tion width of 30.2 MHz. Modelling of the lineshape using 
the optical Bloch equations (OBE) indicates that the spec- 
tral width is limited by the laser linewidths, which are ap- 
proximately 1 MHz. Beam Ai is collimated to a waist of 
1.07 ± 0.02 mm, while beam A2, is focused to a waist of 
12.3 ± 0.2 fim. Hence, a narrow column of Rydberg atoms 
are created in the cold atom cloud, as seen in Fig.[T](b). 

To probe the spatial distribution of the Rydberg atoms, 
the second valence electron is excited using a laser tuned 
near the Sr+5si/2 — > 5p3/2 transition at A 



408 nm 



(Fig. \T\ (a)). For low angular momentum Rydberg states, 
the doubly excited atom ionizes very rapidly (--^ 100 ps). 
The excitation spectrum of the second electron is thus 
broadened to approximately 50 GHz at n ■~ 50 f^J]- We 
drive this transition using an extended cavity diode laser 
stabilized to a wavemeter The autoionizing laser A3 is fo- 
cused to a waist of 10 fim and is pulsed on for 1 fj,s di- 
rectly after the excitation beams have been extinguished to 
prevent direct photoionization. The autoionizing laser is 
orthogonal to Ai and A2. The ions that are produced are 
directed towards a micro-channel plate (MCP) via a weak 



electric field pulse (2.9 V cm~^) that is applied after the 
autoionizing beam is extinguished. The electric field pulse 
is insufficient to field ionize the Rydberg atoms. 

An example of the ion signal generated by autoionizing 
the Rydberg atoms is shown in the inset to Fig.|2](a). Indi- 
vidual ions arriving on the MCP generate a voltage pulse. 
The individual voltage pulses are 2.5 ns wide and are dis- 
tributed over an envelope of 2 fis due to differences in ar- 
rival time and stochastic processes within the MCP. At low 
count rates, individual ion events can therefore be resolved 
and counted using a high-bandwidth digital oscilloscope. 
The ion count at each spatial position is collected for each 
of the 250 repetitions of the experiment. 

The focused autoionizing laser is automatically scanned 
along the y-axis (long axis) of the ensemble by moving 
the focusing lens assembly using a high-precision trans- 
lation stage. Using laser interferometry we measure an 
overall positional accuracy of 2 ^m for a 1 mm displace- 
ment, which is significantly less than the laser spot size. 
By adjusting the angle of the scanning beam, the beam can 
also be moved in the x-direction, enabling us to build up a 
2D-image composed of ID-slices. The focal plane is im- 
aged onto the CCD camera that is used for the fluorescence 
imaging, providing an additional check of the beam dis- 
placement, and a direct calibration of the position of the 
laser focus relative to the ground-state atom cloud. 

A 2D plot of the mean ion signal as a function of position 
is shown in Fig.[T](c). In the absence of background counts, 
the autoionization signal is proportional to the number of 
Rydberg atoms at the intersection of the excitation and au- 
toionizing beams. The image in Fig. [T] (c) therefore rep- 
resents the density distribution of Rydberg atoms in the 
sample. Along the y-axis the Rydberg distribution is de- 
termined by the ground-state spatial distribution, as can 
be seen by the close agreement between the ion signal 
and a slice through the corresponding fluorescence image. 
Along the x-axis a much narrower distribution is measured 
(FWHM = 32 ± 2 ^m), arising from the tightly focused A2 
beam. 

A closer look at the mean ion count for a single slice is 
provided in Fig. |2](a). The small, but non-zero, ion count 
rate in the extremities of the distribution is due to sponta- 
neous ionization of the Rydberg atoms [28]. Spontaneous 
ionization can occur anywhere in the Rydberg ensemble, 
unlike autoionization which only occurs at the intersection 
with the autoionization beam. Even so, at the centre of the 
ensemble the ratio of autoioniziation to spontaneous ion- 
ization is w 20:1, showing the signal-to-noise ratio of the 
technique is very good, and justifying our previous disre- 
gard of the background signal. We have also verified that 
the density of background ions is far below that required 
for an ultracold plasma to form [|llll27ll . 

As well as the mean ion signal, this number-resolved de- 
tection method also provides spatially resolved measure- 
ments of the statistical distribution of counts. Figure. |2](b) 
shows the complete counting distribution for the mean sig- 
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FIG. 2. (Color online) (a) Mean ion signal as a function of y-axis 
position for a single slice. Error bars indicate the standard error. 
Inset: Output signal from the MCP, showing single-ion pulses 
that are counted using a threshold (red line), which excludes noise 
and ringing, (b) The statistical distribution of counts correspond- 
ing to the slice shown, plotted as the sum of the number of repe- 
titions Ni in which i counts were detected. Inset: The histogram 
of Ni at the position indicated by the red line. The black dots with 
error bars show a Poisson distribution with the same mean count 
rate. 



in the wings, the fluctuations are clearly super-Poissonian 
(Q > 0), and the Q-parameter follows the shape of the 
mean ion count. The origin of these fluctuations is techni- 
cal noise in our experiment; a simulation (Fig. |3](b)) that 
includes the independently-measured, shot-to-shot varia- 
tion in atom number, laser detuning and intensity is in good 
agreement with experiment. Based on recent calculations 
of the interaction strength [31], we do not expect to observe 
the Rydberg blockade at the densities that can be reached 
in the relatively hot, high-scattering rate conditions of our 
magneto-optical trap. Indeed, we have seen no evidence of 
excitation suppression or sub-Poissonian statistics for Ry- 
dberg states up to n = 75. 
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FIG. 4. (Color online) (a) The fitted amplitude S and (b) 1/e^ 
width of the Rydberg spatial distribution along the x-axis as a 
function of autoionizing laser power. The red solid lines are the 
results of the OBE simulation described in the text. 
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FIG. 3. (Color online) (a) The Mandel Q-parameter as a func- 
tion of y for the data shown in Fig.|2](a). (b) Simulated Mandel 
Q-parameter, based on the same mean count rate and including 
technical noise. The red lines indicate Gaussian fits to the corre- 
sponding mean ion signal. 



nal slice shown in Fig.|2](a). At each spatial sampling point 
we obtain the histogram of the number of counts, as shown 
in the inset. Access to the statistics, as well as to the mean 
is a major advantage of this technique since the onset of the 
quantum many-body Rydberg blockade regime [4] is asso- 
ciated with sub-Poissonian excitation statistics [29]. From 
the statistical distribution shown in Fig. |2] we have calcu- 
lated the Mandel Q-parameter, which characterizes the de- 
viation from a Poisson distribution llsoll . A plot of the spa- 
tially resolved Q-parameter is shown in Fig. |3](a). Except 



The spatial resolution depends on both the achievable fo- 
cusing of the autoionizing beam, and on its power. At a 
spot size of '-^ 10 ^m, a few hundred microwatts is suffi- 
cient to saturate the autoionization probability for a 1 ^s 
pulse. As the laser power is increased, this saturation leads 
to a broadening of the effective detector size, as autoion- 
ization becomes appreciable even in the wings of the laser 
beam. The saturation of the detection efficiency and the 
corresponding broadening are shown in Fig. 5] Using the 
OBE, and the known transition dipole matrix elements and 
decay rates, we have constructed a numerical model of 
the excitation and detection process that includes the fi- 
nite beam sizes. As shown in Fig. |4] this model is in good 
agreement with the data. There are two fitting parame- 
ters: the overall detection efficiency, which we find to be 
0.21 ± 0.04; and the zero-power width of the Rydberg 
distribution, which we find to be 32 ± 2 ^m. The latter 
is slightly broader than we would expect from the mea- 
sured beam waists of the A2 and A3 lasers - possibly due to 
the laser beam mode quality after transmission through the 
vacuum chamber viewports. 

By exploiting the absence of interaction effects, and the 
large polarizability of Rydberg states, we can apply scan- 
ning autoionization microscopy to the spatially resolved 
measurement of electric fields. Figure |5] shows the effect 
of a spatially inhomogeneous electric field that is applied 
during the Rydberg excitation pulse on the Rydberg dis- 
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FIG. 5. (Color online) (a) Slices through the Rydberg spatial 
distribution at n = 75 with (black circles) and without (blue 
squares) an externally applied inhomogeneous electric field. The 
red solid line is the solution to the simulation described in the 
text, (b) The dashed line shows the electric field extracted from 
the fit to the data above. The solid line shows the Stark shift due 
to the electric field. The colored bars show the standard devia- 
tions of the electric field and Stark shift. 



tribution. The wings of the cloud are Stark shifted out 
of resonance, and the Rydberg distribution and ground- 
state distribution no longer coincide. To extract the elec- 
tric field from these measurements, we use the OBE model 
that fits to the data shown in Fig. [5] The spatial variation 
of the electric field strength is described as a second-order 
polynomial with unknown, adjustable coefficients. From 
this field profile, we calculate the position-dependent Stark 
shift using the polarizability of the 5s75d state ||32ll and 
taking into account the angle between the applied electric 
and magnetic fields. The Rydberg excitation probability as 
a function of y is then calculated using the OBE model, 
with the Stark shift included as a spatially varying detun- 
ing. By multiplying this excitation probability by the cloud 
shape obtained in the absence of electric fields, we gen- 
erate a model Stark-shifted spatial distribution of the Ry- 
dberg excitations. The model distribution is compared to 
the measured data as shown in Fig. |5l and the coefficients 
of the polynomial are adjusted using a least-squares fitting 
routine. An additional fit parameter takes into account a 
possible additional (field-free) detuning. 

For the data shown in Fig. [5] (a), the spatially de- 
pendent Stark shift As and electric field magnitude 
\E\ is shown in Fig. |5] (b). The quadratic term was 
found to be negligible, and we find an electric field 
gradient of 0.78 ± 0.04 V cm'^ with an offset of 
-0.049 ± 0.002 V cm-^ The detuning of A2 from res- 
onance was found to be —2.2 ± 0.2 MHz. The error on 
these measurements includes the error from the fitting pro- 
cedure, as well as an estimated 2.5 % error in the polariz- 



ability. Overall these results are in agreement with a sim- 
ple model of the electric fields produced by the electrodes, 
but the precision of the measurements is significantly better 
than the accuracy of our calculations, which are limited to 
± 20 % due to uncertainties in the geometry and dielectric 
constants of the materials inside the vacuum. 

In conclusion we have utilized the autoionizing process 
in strontium to map the spatial distribution of Rydberg 
atoms in a cold cloud of ground- state atoms. We note that 
although a non-resonant focused laser beam could be used 
to photoionize any Rydberg atom, the ratio of the oscil- 
lator strength for the autoionizing transition relative to di- 
rect photoionization is w 10^^ for the 5s56d state and 
scales as |33], hence using autoionization makes the ex- 
periment tractable. Although this work focused on stron- 
tium the method is applicable to any multi-valence elec- 
tron atom or molecule with autoionizing resonances. The 
method is sensitive enough to detect single Rydberg excita- 
tions, and provides a spatially resolved measurement of the 
excitation statistics. In addition, the efficiency and resolu- 
tion are independent of externally applied electromagnetic 
fields, enabling the high- sensitivity, high-resolution map- 
ping of a weak, spatially inhomogeous electric field. The 
spatial resolution is set by the spot size of the autoionizing 
laser, and could be straightforwardly reduced to 1 fim using 
a higher numerical aperture, which is significantly smaller 
than the correlation length in strongly-interacting Rydberg 
gases [6]. In combination with the single-atom sensitiv- 
ity and number resolving capability that we have demon- 
strated here, scanning autoionization microscopy therefore 
also provides an ideal tool for probing correlations in cold 
Rydberg gases. 
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